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DRAFT 


Recently there has been developed a reformulation of General Relativity 
- referred to as the null surfaee version of CR - where instead of the metric 
field as the basic variable of the theory, families of three-surfaces in a four- 
manifold become basic. From these surfaces themselves, a conformal metric, 
conformal to an Einstein metric, can be constructed. A choice of conformal 
factor turns them into Einstein metrics. The surfaces are then automatically 
characteristic surfaces of this metric. 

In the present paper we explore the linearization of this null surface theory 
and compare it with the standard linear GR. This allows a better understand¬ 
ing of many of the subtle mathematical issues and sheds light on some of the 
obscure points of the null surface theory. It furthermore permits a very sim¬ 
ple solution generating scheme for the linear theory and the beginning of a 
perturbation scheme for the full theory. 


I. INTRODUCTION 

It is our intention to introduce a new point of view to the Einstein theory of gravity. 
We switch our attention from the metric tensor, conventionally viewed as the fundamental 
held of the theory, to a different set of quantities that are to be more basic than the metric 
itself, from which the metric can be extracted as a derived concept. The quantities that are 
our candidates to serve this purpose are null (hyper)surfaces of the spacetime. Null surfaces 
encode the information of a Lorentzian metric up to a conformal factor and the surfaces are 
simply sets of points, that do not carry the elaborate behavior of tensor helds on a manifold, 
therefore being (perhaps) suitable for a more fundamental description of the spacetime from 
which additional structure can be derived, if desired. Our idea is to obtain a description of 
gravity in which we are able to understand Einstein spacetimes in terms of null surfaces. 
In this sense, we depart from the more established view of thinking of GR as a theory for 
the metric as a held tensor Qab on a manifold. We refer to this new approach as a “surface 
theory” of General Relativity. 

The essence of the null surface approach to gravity (hrst introduced in 1983 [|l|) is the 
use of a sphere’s worth of null coordinate systems which allows (up to a conformal factor) 
for all the components of the metric tensor to be obtained by only differential and algebraic 
manipulations on a single function. The idea is to introduce three parameters to label 
null surfaces in the manifold, i.e. a two parameter family of foliations. Since every one- 
parameter family of surfaces can be used to dehne a coordinate system (e.g. the past light 
cones emanating from a world line), then we are actually enlarging our set of variables by 
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adding two more variables to the coordinates. These two additional variables can be taken 
as coordinates on the unit sphere. 

More specihcally, we begin with a Lorentzian manifold M with a metric Qab, a given 
three parameter family of null hypersurfaces on the manifold, namely a function of the form 
F(a;“, u, (, () = 0, the three parameters being u, ( and (, or alternatively hj u = Z{x°', (, () 
with g°‘^Z,a Z,b = 0. Then every set of null surfaces with hxed value of ( and ( (complex 
stereographic coordinates on 8“^), a foliation, can be chosen as the level surfaces of a null 
coordinate u. If we apply an 9 operation (essentially to the function Z, we obtain 
another function of a;“, ( and ( which for hxed ( and ( is in general independent of the 
previous one. We dehne this as a second coordinate u. In the same fashion, the remaining 
two coordinates of a new coordinate system for hxed C and ( can be dehned by applying 
8 and 99 to the original function Z, and will be denoted a) and R respectively. What 
we obtain is a new coordinate system for every value of the parameters ( and (, i.e., a 
sphere’s worth of coordinate systems. The structure we have is that of the bundle of null 
directions over M, with each hber an and a coordinate system associated with each hber 
point. There is a major advantage to having this “bundle” of coordinate systems: one can 
diherentiate the components of any tensor with respect to the parameters, i.e. vertically, 
and obtain relationships among the components themselves. In other words, the components 
of a tensor will not be independent of each other, and it will be feasible to obtain most of 
them by knowledge of only one for all values of (C;C)- This and the fact that one of the 
coordinates is null provide the basis for describing the metric tensor completely in terms of 
the single function Z(x“, C, (), which becomes a basic variable of null surface theory. 

The dynamics of the function Z has been a troublesome issue ever since its appearance. 
In principle, all 10 of Einstein equations are encoded in only one equation for the con¬ 
formal factor n and Z{x°‘, () by virtue of the sphere’s worth of coordinate systems. There 

are, in addition, several “auxiliary” conditions that guarantee that the metrics, eventually 
extracted from Z{x°‘X)C)) cire all diffeomorphic to a metric gab{x^)- Considerable success 
was originally attained by introducing the holonomy operator around closed loops as an 
auxiliary variable |^, though the resulting theory was quite complicated. Recently a great 
simplihcation was achieved where the full GR was expressed as equations for Z and hi. 
There however still remained certain difficult conceptual issues, which led us to study the 
linearized version developed below. This gives a new insight into the structure of the full 
equations. 

We hrst review, in Sec.II, the conventional approach (conventional as opposed to the 
null surface approach) to linearized gravity, showing the Einstein equations in terms of the 
hrst order correction to an assumed underlying hat space. We make a distinction between 
trace and trace free part of the linearized metric and give the equations in terms of these 
two quantities instead of the full hrst order tensor. We then show how a choice of gauge, 
in the conventional sense, can lead to a new coordinate system that is similar to the one 
described above, except for the critical (C, () dependence. We believe that this comparison 
of the conventional view with the null surface view clarihes many of the issues of null surface 
theory in hrst order. 

A complete discussion of hrst order null surface theory for gravity with sources is found in 
sections III and IV. The attention is focused on the role of sources and a possible cosmological 
constant on the two relevant quantities. We hnd that all the components of the metric can be 
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expressed algebraically in terms of derivatives of a single spin-weight-2 function denoted A, 
which is dehned by A = and the trace of the metric v which to hrst order is equivalent 
to a small correction to the conformal factor It is remarkable that the metric does not 
depend directly on Z but only through A, which means that the metric does not know of 
the hypersurfaces themselves, but of a related quantity A whose geometrical meaning will 
be discussed elsewhere. 

II. THE LINEARIZED FIELD EQUATIONS 

We begin with a spacetime and a metric Qab that differs from a flat metric by a small term 
with Qah = Tjah — lab- & perturbatiou expansion to hrst order, the Einstein equations 
for Qab yield a system of linear second order equations for ■jab on the hat background rjab- In 
particular, the linearized Ricci tensor becomes @]: 

Rab = -d''d^bja)c +^D'Jab + ^dadbj ( 1 ) 

R = + Dq (2) 

where da is the hat derivative, □ = d^^dc, R and 7 are the (linearized) traces of the Ricci 
tensor and the hrst order correction to the hat metric. Indices are raised and lowered with rjab- 
We decompose the hrst order metric into its trace-free part and its trace: jab = \lTlab + Qab 
and redehne the trace 7 in terms of a scalar u in the following way: 7 = Su. The Ricci 
tensor, in terms of the redehned “trace” u and the trace-free metric Qab, becomes 

Rab = 2dadbiy + r]abOu- d'^d[bqa)c + ^ □ ( 3 ) 

and the Ricci scalar becomes: 

R = 6 □ z/ - d^dUad. (4) 

We assume that the source of the Einstein equations Rab — \gabR = 2Tab is a hrst order 
stress-energy tensor that is conserved, and does not depend on the hrst order correction to 
the metric, but on the hat metric at most. As dehned here. Tab absorbs the factor k/2 that 
is usually explicit. Tab can be written as Tab = tab + \'UabT with tab trace free and T its trace. 
The conservation law d°‘Tab = 0 for the stress energy tensor in terms of tab and T takes the 
form: 


dbT = -4dRab. ( 5 ) 

We will refer to Qab loosely as the metric and to v as the trace, to avoid unnecessary 
complexity in the terminology which may obscure the content of the work. 

The complete linearized Einstein equations in terms of v and Qab read: 
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dadbiy -T]abniy- -d’^d^bQaj + -T\qab + ^Vabd^'d'^qad = Tab- ( 6 ) 
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It is possible to choose coordinate systems so that the flat metric has only the following 
nonvanishing components: r^oo = 2, r^oi = 1, and = —1. A way to do this is by choosing 
a null tetrad (/a,hta) where la and are real, nia = fla + iba with Oq and ba real 
spacelike vectors, and all the scalar products being zero except I’^Ua = —m'^rha = 1. Then 
we would take la, ma, fha and the combination Ua — la as our coordinate basis vectors. 
The inverse has the following non-vanishing elements: = —2, = 1 and = 

— 1. We will understand the equations (|^) as equations for the two quantities u and qab 
in these particular flat coordinates, so that both of them are functions of the coordinates 
,x~) and partial derivatives are taken with respect to these coordinates. An 
example of these coordinates would he x^ = z — t, x^ = z, x~^ = x + iy and x~ = x — iy. 

The linearized equations (|^) admit gauge freedom for both the metric and the trace in 
the following form: 


Qab Qab T Qa^b T ^b^a ^Qab^ 


( 7 ) 


( 8 ) 

where ^a is a free vector held. 

Using this gauge freedom we can make q[f^ = 0 for 6 = 1, — and goi+ Q+- — 0- Notice 

that since q'ai, is trace free then 0 = rj^^^q'ab = ‘^{Qoi~ Qii~ Q+-) implies = q'oi = 0 as well. 
There are no restrictions on the trace z/. These conditions on the gauge can be expressed as: 

di^b + db^i - ^Vib&'^c + qib = 0, 

which can be solved explicitly with a given qib- There is still the remaining gauge freedom 
of the homogeneous system: 

di^b + db^i - ^Vibd^'^c = 0 . 

The solutions to this system leave four unspecified functions of the coordinates x'^ and 
x~, while the dependence on x^ is at most quadratic. The corresponding divergence 
is linear in x^ with two free functions of the remaining coordinates as the coefficients. We 
will return to the issue of the gauge freedom in Sec.IV when we introduce a special family 
of coordinate systems and restrict to the transformations within the family. 

Summarizing, the linearized metric has the form: 


Qab ^ {.Vab Qab') (1 2z/) 


/2 1 0 0 \ 
10 0 0 
0 0 0 -1 
VO 0 -1 0 / 


with inverse: 
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If we now look at the a = 1, 6 = 1 component of the Einstein equations (|^) we obtain a 
single equation that relates the sources to only the trace of the hrst order metric, since all 
the terms that involve the trace free part qab are identically zero by virtue of our choice of 
gauge. The equation will be: 

= tu. (9) 

This equation tells us that the source directly drives u, i.e. u is given in terms of a 
double integral of the source tn plus any solution of the homogeneous equation. Note that 
the trace free component tu is equal to the full Tu- 

The homogeneous solutions to (j^) are all functions that are linear in the coordinate 
with two free functions of the remaining coordinates as the coefficients, i.e. uh = a + 

We will return to this later. 

The motivation for the discussion up to now has been to establish a basis for the under¬ 
standing of first order null surface theory from a more conventional point of view. We 
will see that, in a sense to be discussed, Eq.(^ is the single equation that encodes the 
dynamics of first order Einstein fields from the null surface point of view. The 
new ingredient in null surface theory that explains why the remaining nine equations do not 
give any other relevant information is that it makes exhaustive use of the introduction, not 
of a single coordinate system with these properties, but a smooth two-parameter family of 
them. We will devote the following sections to a rather complete discussion of the ideas of 
null surface theory to hrst order. 


III. A SPHERE’S WORTH OF COORDINATE SYSTEMS 

We begin with a manifold M with a Lorentzian metric Qab (we might also want a boundary 
to M if we wish to discuss asymptotically hat spaces). A choice of coordinates based on 
null surfaces is made in the following fashion. We choose a family of null surfaces that are 
labeled by points of a three parameter space: Rx 3“^, with a real parameter u on i?, and two 
stereographic coordinates (C, C) This is expressed by F(x“, u, C, Q = 0, or preferably 

by M = Z(x“, Cj 0- An immediate meaning to the function Z from this point of view is that 
for every hxed value of the parameters (C, C) the level surfaces oiu = Z are null surfaces i.e. 
g°‘^daZdbZ = 0 , labeled by the parameter u. 

For the construction of our coordinate system we choose four functions of the points 
that are simply Z, dZ, 8 Z, and 8 dZ at an arbitrary but hxed value of the parameters 
(C) 0 - We choose to work with the S and 8 operators instead of ordinary partial derivatives 
because of their covariant nature 0. We thus dehne an S'^’s set of coordinate systems for the 
spacetime in the following way: 

U = Z(x“,C,C) 
o; = 8Z(x“,C,C) 

(I; = 8Z(x“,C,C) 

R = mZ{x^XX) ( 10 ) 
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We assume familiarity with the 6 and 9 operators. See reference [^. 
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where the parameters {(, () have hxed values (after differentiation). As a matter of notation 
we refer to them as 0* = (u, R, uj, uj), with i = 0,1, +, —• We thus have 6** = ((). 

[As an example of this construction in an asymptotically flat spacetime we can think of 
the interior spacetime points as being while the boundary is a three-dimensional space 
(the null surface J'~^) which can be given the coordinates u, ( and (. From a hxed value 
of (uXX), the past light cone is the set of points that satisfy Z[x'^XX) = thus u 
labels past light cones from hxed null generators of . Early history of the null surface 
theory gave Z the name light cone cut function ||l||; we will sometimes refer to Z as the cut 
function. On the given past light cone, the values of (a;, a)) label the null generators, thus a 
particular value of {u,uj,iD) is a null geodesic emanating from the point (n, C, C) The 

last coordinate i? is a parameter along the null geodesics (in general not affine).] 

Two diherent values of {(, () give two diherent coordinate systems. However, whatever 
coordinate conditions are implied by (|^) will be true for all values of (C,C)- Since, by 
assumption, the dependence on ((, () is smooth, one can go back and forth from one coor¬ 
dinate system to a neighboring one by diherentiating with respect to the parameters ((, (), 
or equivalently, by applying 9 and 8 operations. 

From (pffiD we obtain a sphere’s worth of coordinate tetrads. The gradients give us a 
coordinate set of covectors: 


c= 

= 8^,a 
= 8^,a 

= 88Za. (11) 

The coordinate tetrad vectors 6*f dual to 0^^ are dehned by = Sj. 

We will be concerned with the linearized version of the Einstein equations, in which all 
terms are flat second derivatives of the hrst order correction to the flat metric. The function 
Z to hrst order consists of a zeroth order Zm and a hrst order term ; 2 , i.e. 

Z(a:“, C, 0 = Zm{x% C, 0 + ^(^“, C, 0 (12) 

where Zm is supposed to be known and z becomes the hrst order variable. The role of Zm 
is to hx the zeroth order hat background: the hat intrinsic coordinates and tetrad. It will 
not be necessary to keep hrst order terms in the tetrad, since, in what follows, they will 
always appear contracted with tensors that are already hrst order (see remark below). We 
will work with the hat tetrad obtained from the function Zm for hat space dehned by a 
sphere’s worth of null planes, which is easily seen to be [|l|: 

Zm = xRa{C,0 (13) 

where are standard Minkowskian coordinates and 

^“(C,C) = ^|^^^^((1 + CC), C + C, ^(C-C), -1 + Cc) (14) 

is a null vector for all (C,C) with respect to the hat rjab- The index is raised with rjab with 

signature (H-). It is customary to dehne a null tetrad out of this function where the 

remaining vectors are: 
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rUa = dla 

Ua = BS/a + la (15) 

SO that n°‘la = 1, rrf'fha = —1 and the remaining scalar prodncts are zero. Following onr 
approach, we dehne 0* = {u,R,uj,iD) as in (PD]). Then: 

u = xRa{.CX) 

UJ = a;“ma(C,C) 

(h = x'^TflaiCX) 

R = x^MC,C)-UC,C)). (16) 

Onr tetrad is then: 

(^aj ) ( 1 ^^) 

and dnal: 

= (n“ + -m“, -m“ ) (18) 

and the partial derivatives with respect to the coordinates are di = O^da- 

The flat metric rjij in these coordinates has three non-vanishing elements: r^oo = 2, r^oi = 1; 
and = —1 with the inverse having the elements: = —2, = 1 and ri~^~ = — 1 

non-vanishing, i.e., is of the form we nsed in Sec.II. 

We end this section by pointing ont that this (C, C)-dependent set of coordinate systems 
has the properties (i.e. satisfies the same set of coordinate conditions as in Sec.II.) reqnired 
to obtain the eqnation for the trace u of the first order metric from the Einstein system of 
eqnations, namely Eq.(^), dfu = tn, for each valne of (C^O- The advantage now is that 
we have a smooth two-parameter family of snch coordinate systems, so we are allowed to 
perform differentiations with respect to {(, () withont changing the coordinate conditions. 
These remarks will be amplified npon in the next section. 


IV. FIRST ORDER NULL SURFACE THEORY 
A. The null theory variables n and A and kinematic equations 

The components of any tensor referred to the new coordinates are obtained by contraction 
of the old tensor with the basis or dnal tetrad. In particnlar, the metric components with 
respect to the new coordinates denoted by gij are given by gij = gabOfd'j, with inverse 
g^j = g'^%9^^,. 

There is a distinction between ij and ab indices. Qnantities like gab or Tab are “local” in 
the sense that they depend only on the local coordinates of the spacetime x°'. However g^j 
or Tij depend on 0® and (C, C) as well. Similarly, the partial derivatives da commnte with 9 
and B, whereas di do not in general. The important commntation relations of 9 and B with 
di are given later. 
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In these coordinates [Q, all the components of the metric can be expressed in terms of 
two functions. One of these two functions is defined as follows: 


A = 52 ^. (19) 

All the components of the metric get expressed as times a simple algebraic function of A 
and its derivatives. Therefore can be taken as a conformal factor = g^^ = g'^^9\9\ = 
g°'^Z,a Thus the metric is conformally a function of A, the conformal factor being 

g^^. This follows essentially from the fact that g^^ = g°'^Z,a Z,i, is zero by assumption, since 
u = Z = constant is a null surface, and the remainder of the components can be computed 
by applying S and 0 to g^^ = 0 an appropriate number of times, using the non-trivial fact 
that = Bg°'^ = 0. For example, from 

g^^ = g-^Z,aZ,,= 0 ( 20 ) 

and 

g^+ = g^^Z,adZ,,, (21) 

by applying 0 to g^^ = 0 we obtain 

^ 2 ^ 0 + ^22) 

and therefore g^~^ = 0. Applying dQ{g‘^^Z,a Z,h) = 0 we get 005 f°° = 2g^^ + 2g^~ = 0 and 
therefore 

9^- = (23) 

Taking Q'^g^^ = 0 gives the expression of g in terms of g^^ and A. Taking 00^ of g^^ = 0 
gives g^^. Finally, by taking 0^0^ of g^^ we obtain g^^. The complex conjugates of all these 
operations yield the remaining components of the metric. 

The linearization can be carried out by noticing that Zm = x°'la satishes 

Q'^Zm = Q^Zm = 0 . 


Thus the zeroth order does not contribute to A and therefore A is a hrst order quantity. By 
dividing the g^^ by g^^ and linearizing in A we have g'^^/g^^ = + /*-^ where is first order 

and has the following hve non-vanishing elements: 

/II = -^02Ai + 0A+ = -^0'Ai + 0A- 
/i+ = _l8A^ 
f- = -^SAi 

/++ = -Ai 

/-- = -Ai (24) 

The non-vanishing fij are the following: 
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/oo = + SA+ = -^9'Ai + 9A_ 

/o+ = ^BAi 

/o- = ^BAi 
/++ = -Ai 

/__ = -Ai (25) 

where we have adopted the notation c?jA = Aj. Note that the hrst order trace-free metric is 
given in terms of A and its derivatives. 

In deriving the components of the metric from = 0 we used repeatedly = ^g°'^ = 
0. This is the statement that all the (C, C)-dependent metrics obtained by g'’^ = g°‘^9'‘^9\ are 
in fact diffeomorphic to g°'^, i.e., they are all equivalent. The conditions 'dg°'^ = Qg°'’^ = 0 
constitute a set of 20 equations of which some, when expressed in these coordinates, give the 
metric (Eq.(|25|)). Some are identities. The remaining equations, given below, are referred 
to as the “metricity conditions”. We discuss their meaning and signihcance in [^], including 
an explicit derivation of the metricity conditions from a slightly different point of view. 

Notice that / = fij = 0 automatically. On the other hand, the conformal factor 
g^^ = 1 + 2v is equal to = 1 to zeroth order plus a hrst order correction 2z/ which 
can be identihed with the trace of the full hrst order departure from hat space, since g''^ = 
^01 _j_ jij'^ _ ^ij _j_ _|_ jv ancj therefore the trace of the hrst order correction is 

indistinguishable from the hrst order conformal factor. In fact, we will make no distinction 
between the trace u of section II and the hrst order conformal factor of the current section. 

We can continue to take 9 and 8 of the components of the metric that we have already 
obtained. The only new relations (the metricity conditions) that we can obtain in this 
way P are the following: 

(*) By writing g^^ = g°''^Z,aQQZ,h, applying 9 to g^^ and linearizing in A we obtain the 
relationship between u and A: 

29z/ = A+ + ^BAi. (26) 

(A) If we apply 9^ to g^^ = g°'^Z,a Z,!, = 0 (the 6'+6'’{)-component of Qg^"^ = 0) and linearize 
in A we obtain the following relation for A only: 

0 = A_ - ^9Ai. (27) 

(B) We can now take QQ^g^^ = dd^{g°'^ZaZh) = 0 (the ^(^^■{(-component of Qg'^^ = 0). 
After linearizing in A we obtain: 

0 = B(9Ai-2A_). (28) 

(C) Finally, by taking = 0 (the 6')j6'),-component of 9(7“^) we obtain another 

linearized condition on A: 

0 = d^Ai - 29BA+ + 49Ao - 4A+ - 2 BAi (29) 
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(30) 


which, with the help of (p6D gives: 


0 = + 4gAo - SQu - 


These four hrst order relations (*), (A), (B) and (C), are of a kinematical nature in 
the sense that they hold for any linearized Lorentzian metric independent of the Einstein 
equations. They are very important to the null surface approach, and as such they deserve 
a name of their own: we refer to them as the metricity conditions. If the null surfaces 
u = Z(x“, C, C) of a given metric are known and A is computed, then the metricity conditions 
are identities. But if the null surfaces are yet to be found, then the metricity conditions 
become the equations that yield A, from which the metric and, if desired, the null surfaces as 
well, can be computed (the metricity conditions are the requirement that the sphere’s worth 
of metrics be really equivalent). It should be noted that in hrst order theory Eq.(B) 
follows from Eq.(A) i; this leaves only three relevant metricity equations (*), (A) and (C). 
They play an important role later. 

Before we further analyze these kinematic relations and properties of the we need the 
commutation relations between 9 and the di. Since WiF = d{9fdaF) = dO^daF + OfdaQF 
and since F will be chosen as hrst order then: 


5-5 = dd- 

d+d = + do — 2di 

d^d = ddi + d_ (31) 

with the complex conjugate relations between 8 and the dt. In obtaining these relations we 
used the following: 

80° = gz = 0+ 


80+ = g^z = 0 


80 - = ggz = 0 ^ 


80^ = g2gz = -20+. 


(32) 


or alternatively relations ([T5|), ([T7| ) and (pISl). 

We mention here some very useful results that we will need later, namely the divergence 
d^ fij and the double divergence d^d^ fij. With the use of the commutation relations and the 
kinematic relation (*) we can evaluate d^fij. First we notice that since fij = 0 then the 
contraction d^ fij has only three components. Also, since d^ f+j is the complex conjugate of 
d^f-j we only need to compute the components z = 0 and i = —. We display below the 
results and leave the detailed calculation to an appendix: 


remains an identity even in the full theory (see Q) 
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d^f-j = 2didu 


(33) 


with the complex conjugate: 


(34) 

and hnally 

d^fQ. = -2dimv + 2d+^v + 2a_Bz/. (35) 

The double divergence fij turns out to be simply: 

d^d^f.. = -2dlBdu. (36) 

Note that although fij are functions of A the divergences are functions of only u. 

B. The dynamics of u and A 

So far the description of null surface theory has been completely kinematical. Our 
variables (n and A) must satisfy the metricity conditions (*), (A) and (C) to dehne for us a 
Lorentzian metric. We now address the problem of hnding a metric that in addition satishes 
the Einstein equations. The Einstein equations are implemented in the following way. We 
contract the linearized Einstein equations (|^) with (= 6i6\): 

ri^(dadbi' -r]abOu- ^d''d^bqa)c + ^Oqab + ^Vabd^'d^qcd) = Tabl^'l'". (37) 

Because our gauge is the one described in section II, and because is null, we immediately 
obtain: 


ri^dadbP = Tabn^ 


or 


c)P(9‘,C,0 = iii(etC,C). 


(38) 


This is the dynamical equation for A and v. Though it resembles Eq.(^), its meaning is 
different. Because of the (C,C) dependence of it is not a single equation, but an S'^’s 
worth of equations instead. It contains all the equations that make up the trace-free part of 
the Einstein system. 

[Another way to see that Eq.(^^ encodes the trace-free part of the Einstein system is 
the following. We can apply 5 and 8 to Eq. (|38D an appropriate number of times so that we 
obtain the remaining eight of the trace-free equations. It is clear that we can do so because 
til = Tabl°‘l^, and thus an 8 operation on tu goes through Tab and is applied to directly, 
giving a different tetrad vector (m“), and therefore giving a different component of Tab- See 
the appendix for a more complete discussion of this issue.] 

It may appear that we are lacking a tenth equation, the trace of the Einstein system. It 
is true that the trace is lost when we contract the Einstein tensor with However, the 
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information contained in the trace of the Einstein system is redundant (see Eq.(^)). If we 
put the Einstein equations in terms of the trace and trace-free part of the Einstein tensor 
Gab = G^b + \ 9 abG they become: 


G_ab ^ab 


and 


G = T. 


In addition, the Bianchi identities \/°‘Gab = 0 become: 

VfeG = -4 v“ G^b- 

Therefore the trace of the Einstein tensor is determined (up to a constant) by the Bianchi 
identities. However, the Bianchi identities are of a kinematical nature, in that they are true 
for any Lorentzian metric. This means that the trace of the Einstein system is redundant 
when and G are expressed as functions of a metric. Now returning to the null surface 
point of view, we claim that any solutions v = to Eq. (|55D and the metricity 

conditions as well (since the metricity conditions guarantee a metric), will provide us with 
a Lorentzian metric that satishes the complete set of Einstein equations. It is in this sense 
that we say that Eq.( 


encodes the full Einstein equations. 

[We mention that the only piece of information that we lack in the null surface approach 
is the integration constant of the Bianchi identities. This integration constant is the cosmo¬ 
logical constant in the vacuum case, or part of the stress tensor in the case of non-vanishing 
sources.] 

In the remainder of this section we will clarify the issues just raised. First, we will return 
to the conventional form of the Einstein equations and calculate the trace of the Einstein 
Eqs.(|^). Then we will show that in the null surface approach, this trace equation is satished 
if Eq.(PBD and the metricity conditions are satished. Thus, we will have proved our claim 
that we do not need the trace equation. In the null surface approach, the tenth equation is 


thus encoded into the metricity conditions and (|^) . 
Contracting Eq.(H) with 7]°"^ we obtain: 


-3Du+-ffd^f,j = T. (39) 

By virtue of our previous calculation for the double divergence of the metric fij, i.e. Eq.(P^), 
Eq.(^^ turns out to be an equation for only the conformal factor z/, decoupled from A: 


-3 □ z/ - = T (40) 

where T is uniquely given by daT = —Ad^tab up to a constant. This would be the tenth 


Einstein equation. We now show that (^) is satished by virtue of (^). With some work 
(by essentially using the commutation relations (PT|)) it can be shown that Eq.(|55D implies 
the following: 


di{-3Du -d^mu) =diT 


(41) 
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(where diT is equal to —AdHij). Furthermore, we show in the appendix that the metricity 
conditions (*) and (A) imply: 


5(-3 □ z/ - = 0 (42) 

8(-3 □ 1/ - dlmu) = 0. (43) 

These are sufficient conditions to guarantee that c?i(—3 □ z/ — S^BSz/) = diT for all i because 
the application of 8 and 8 on di{—3 □ z/ — dfddu) = diT yields the remaining equations 

if (—3 □ z/ — dfBdu) is a function only of as implied by (|4^) and (|43|) . Therefore any 

solution to Eq.(p8D and the metricity conditions satishes the tenth equation automatically. 
Note that (|4^) and (|4^ imply severe conditions on a and [3 in the solution to (|3^) (see 
Sec.IV.D). 


C. The solution to the metricity conditions 


Our point of view is now as follows: without yet worrying about the details (see discussion 
at the end of this Sec.) the conformal factor z/ can be found explicitly as a double integral 
in R of the given stress tensor fn. Once v has been solved for, it enters the equations for A 
(Eq.®, (0) and (pOD) as a known source only. Our task is then to solve the three kinematic 
equations for A. [It should be noted, however, that the decoupling of the equations for z/ 
and A can not be carried out in this way in the full theory, in which the metricity conditions 
and the dynamical equation are strongly coupled 0.] 

We begin by manipulating (*) and (A) to obtain two very important relations for A. The 
following is merely the derivation of Eqs.( 


and (^61) which we include at this point for 
the sake of completeness (too frequently in this work the results are not easy to reproduce, 
due in part to the non trivial commutation relations ([31|)). 

By applying 8 to (^Tf) and commuting 88 on the right side we obtain: 


8A_ = -88Ai + 2Ai 

By applying now di to both sides and commuting 9i8 in the left side we obtain: 


Wik. + a+A_ = -ai88Ai + 2diki. 


If we commute now 89_ in the left side we obtain: 


9_8Ai — c^qAi + 29iAi + (?+A_ = ^c}i88Ai + 2diKi. 

The c?iAi cancel out from both sides. We now use the basic equation (p^) to eliminate 
8Ai in the left side and obtain: 


a+A_ 


-aiAo--ai88Ai+4a_8z/. 


(44) 
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Next we carry out a similar procedure on First we apply 5 and then to both 

sides. By commuting 5i9 on the left side obtain: 

S9iA+ + 5_A+ = —^9 iS8Ai + 

Further, by commuting Qd+ on the left side we obtain: 


9+9Ai — (?oAi + 2diAi + c}_A+ = —-9i89Ai + 2diQ^u. 


Using 


to eliminate 8A_ on the left side we are left with: 

39_A+ — c^oAi + 29iAi = — ^9 i58Ai + 2diQ^h'. 


Using d+A- = d-A^ and inserting d+A- from we obtain our hrst basic result : 

49iAo + 9 i88Ai — 2diAi = 129_8z/ — 29i8^z/ (45) 

Eliminating c}i88Ai from via (]4^) we obtain our second basic result : 

□A = 2{dodi — didi — d+d-)A = 45_8z/ — 2did‘^u. (46) 


A thus satishes the wave equation with the conformal factor playing the role of the source. 
[This equation, (|46|), has been derived here just for completeness and will not be used. Its 
use will be presented in a future paper where A will be the basic variable and Z purely 
auxiliary. In the present paper we use A as auxiliary while Z plays the primary role.] 

Equation (^5]) can be integrated on R, assuming asymptotically flat boundary conditions, 
to yield the following: 

Ao = - 78 §Ai + ^Ai + 3 /" 9_8z/ci.R-^8V + a,o(M,C,C) (47) 

4 Z J oo Z 


where a{u, (, () is free radiation data. 


Taking 8 of Aq and inserting it into the third metricity condition (|30D we End: 


rR 


- 8 '^Ai = 8cr,o+38 / d-'dvdR — 2^v - 8 '^ 8 i/ 


(48) 


and its complex conjugate. 

By applying 8 to (|30D we obtain: 

U 
-( 

4 

Now using (^8|) and its complex conjugate in the previous equation leads to: 


82An = --m^Ai - -m^Ai + 88z/ + 8888i/. 
4 


rR 


rR 


82Ao = 8V,o+82^,0+382/ d_QjydR + 3d^ a+8z/di? - 388i/- 288^8// 


(49) 


Integrating once in u and using A = hnally yields the fundamental equation for the 
cut function Z: 
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+ g 2 - ^ + Sd^J^d+dudR - 3ddu - 2m^Qiy^ du. (50) 

Equation (|^) , our final equation, equivalent to the linearized Einstein equations, can be 
solved by means of a simple Green’s function for the operator 0^9^ (essentially the double 
laplacian on the sphere). The Green’s function is given explicitly in appendix G |^. As 
complicated as it appears, the right hand side is simply the source, a known function of 
(a:“, C, Q once the conformal factor has been solved for. 

To solve for the conformal factor in terms of only the source tn is a task more subtle 
than it appears at first sight. Formally, one can write the solution to Eq.(^8D as the double 
integral of tn plus a homogeneous solution that is linear in R and that contains the values 
of V and c?iz/ at a fixed value R = Rq, as shown earlier. The homogeneous part is severely 
restricted, however, by the trace equation. The trace equation, in terms of a and fd becomes 
quite cumbersome in the case of general (arbitrary) sources, involving the values of tn and 
different derivatives of tn in a non-trivial manner at the surface R = Rq. We will bypass the 
problem of solving for the homogeneous part of the conformal factor in the general case by 
restricting ourselves to the case where there exists a surface R = Rq on which tn vanishes 
as a function of {u, Ro,uj,ujXX)- This includes the asymptotically fiat case (where tu = 0 
at J'~^) but is slightly more general. 

For the particular case of sources that vanish at a surface R = Rq, the conformal factor 
can be solved for unambiguously by: 


u = 



rR 


tiidR'dR 


' Ro 


( 61 ) 


in the gauge such that a = j3 = 0 (see next section for a discussion of gauge). 

The corresponding equation for the cut function is the same as (^) with the lower limit 
of the integrals being replaced by Rq. 


D. Gauge conditions and remaining freedom 

This section is supplementary to the linearized null surface approach described in the 
previous sections. We justify particular choices of the conformal factor u made in past 
works (and in (|^ ) and provide some insights into the issue of the gauge freedom. 

We mentioned before. Sec.II, that by imposing qu = 0 we were still left with gauge 
freedom of the form: 

dXi + dXi - = 0 . 


For i = 1 we have: 


■Ciu — 0 


therefore .^i = (, (). For i = + we have: 


■^1)+ + 'C+u — 0 
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therefore = —0,+ -R + r+ where r+ does not depend on R. Similarly for i = — we obtain 
= —0,_ R + t_ where r_ does not depend on R. The component can be solved for in a 
similar fashion from ~ = 0 but it will not be needed in what follows. The 

divergence is equal to —2(^+,_ +^_,+ ) and therefore = 40, +_ R — 2(r+,_ +t_,+ ). 
At this stage it is clear that the gauge freedom in v (see (|])) is linear in R with two free 
functions of the remaining coordinates (and (C^C)) ^-s the coefficients. 

Additional gauge conditions arise from our (^, 0)-dependent coordinate systems; see 
Eqs.([T0|). The held is just an inhnitesimal variation of the coordinates to neighboring 
coordinates . Explicitly: 

v^ = uRbu = uR^ 
cn' = cn + hen = a; + 
uj' = ijj R buj = u) \ 

R1 = RRbR = RRi^. ( 52 ) 

Since 0'* and for i = 1, +, — are obtained from Q'^ and (P by the application of 9, 8 and 
88, the components and can also be obtained by applying 8 and 8 to ^^ i.e.: 

r = 

= 88^°. 

Thus, r+, r_ and 0 are not independent of each other. If we lower the index with 
we hnd the following: 


6 = ^° = 0 (m,o;,o;,C,C) 

= -r = -8^° = -86 
= -8e° = -86- 

Since = 0 (m, a;,(h, C,C) and the coordinates u, cn and eh depend on (C, C) then taking 8 
corresponds to taking a total derivative with respect to (C^C)- Assuming that 0 is a hrst 
order quantity we have: 

8^1 = 80 = 0,0 80° + 0,+ 80+ + 0,_ 80- + 8'0 
= 0,0 en + 0,_ R + 8^0 

8^^ = 80 = 0,0 800 + 0, +-80+ + 0,_ 80- + 8'0 
= 0,0 w + 0,+ R + 8^0 

where 8' and 8' are taken holding 0* hxed and the dehnitions of en, eh and R have been used. 
Therefore 


r+ = -0,0 et) - 8'0 
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r_ = —0,0 0; — 

With these expressions the divergence is given as: 

d^^k — 40,H— R + 40,0 +20,0— o) + 20,0+ to + 20^0,_ +29^0,+ . (53) 

Here 0, an arbitrary function of {u,uj,lD, (, (), represents our full gauge freedom. The diver¬ 
gence d^^k can be used to modify u, as noted in Sec.II. Recall (^: 

By virtue of Eq. (|38D , u is determined by tn up to a linear function of R. We are interested 
in seeing whether we can gauge away the linear (homogeneous) part of u. We claim that if 
u satishes Eq. (pSD and the metricity conditions, then we have just enough freedom to gauge 
away the homogeneous part of u, in the case where the cosmological constant vanishes and 
there is a surface R = Rq on which tn vanishes. In particular, u would be pure gauge for 
the vacuum case with vanishing cosmological constant. This would justify the choice z/ = 0 
made in earlier works on the linearized vacuum null surface theory [Q. 

The proof of our claim is in the following argument. As was previously shown, any u 
that satishes both Eq.(p8D and the metricity conditions will consequently satisfy 

-3 □ z/ - dlSBu = T, 


the trace equation. Although not needed in general, it proves particularly useful in this 
section. If there exists a surface R = Rq in which tn vanishes, when we substitute u = 
Iro jRQtndR'dR + a + j3R (the general solution to Eq.(^)) into the trace equation, the non- 
homogeneous part of v drops out and we obtain an equation for only a and (3 that has the 
following form: 

80 - 4a,+_ +20,_ u) + 20,+ oj + 28'0,_ +25'0,+ = -A (54) 

where 9' and 9' are taken by holding the coordinates hxed, 0 = 0,o and A is the cosmo¬ 
logical constant. 

One technical point that was needed in the derivation of (|54D is the commutator of 9 
through the integral sign, given by (boundary terms vanish for our particular restriction on 
Hi): 


cR 


i-R 


RrR 


9 / FdR = / dFdR + 


' Ro 


' Ro 


d_FdR!. 


' R& Rq 


In this sense, we think of the trace equation as an equation for just the homogeneous 
part of v. It is remarkable that d^^k (given above (l53|) of the form a + 0i?) is of the form so 
that a and 0 solve the trace equation. In other words, the expression for d’^^k can be chosen 
as the solution to the trace equation (with vanishing cosmological constant). Therefore, we 
can always choose a gauge so that the conformal factor u has no homogenous part, if the 
cosmological constant is zero and the source Hi vanishes on a surface R = Rq. 

The solution to the trace equation with non-vanishing A, modulo gauge, is 


u = \{uR — R^ — ujuj). (55) 

This solution can not be gauged away. We are free to choose A to vanish or not to vanish. A 
general way to rule out a possible cosmological constant is by requiring asymptotically flat 
boundary conditions for the conformal factor. 
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V. CONCLUDING REMARKS 


The null surface theory approach to GR to hrst order just described is meant to serve 
as a reliable test of the concepts of the full null surface theory. We found it extremely 
useful in clarifying the full version of the theory, which will be published elsewhere 0. As a 
restriction to the linearized regime, however, it’s utility is limited and some results can not 
be carried over to the full theory. For example, we found that condition (B) in Sec.(IV.A.) 
is an identity in first order, but the fact that it remains an identity to all orders can not be 
established from the linearization. Nevertheless, it gave a hint as to what to expect from 
the full theory with respect to the relevance of (B). 

The linearized null surface theory departs conceptually from the full theory in that it 
is a theory for the quantity A on a fixed known background. In this sense, many of the 
complications of the full theory do not have a place here, in particular, the difficulty that 
the coordinates and A are conceptually intertwined in the full theory, which makes the 
equations extremely lengthy (and the ideas at the start quite confusing). It’s a fact that 
the linearization of the null surface theory becomes itself a separate theory, easier to handle, 
where the methods used do not need to agree completely with the ones needed in the full 
theory. For example, the conformal factor p can actually be solved for without knowing A, 
and then it can be used as a known source to find A, a fact that does not apply to the full 
theory. 

Throughout this work, the null surfaces are known to zeroth order and get corrected 
by the hrst order A, i.e., A gives a hrst order correction to the null planes introduced in 
Sec.II. One can develop a perturbative scheme in which the surfaces are null planes in hat 
space, and get corrected successively with every higher order step in the hierarchy. The 
complications in the equations increase considerably at the second order level, making it, at 
least for the moment, a very difficult task to accomplish explicitly . A formal procedure can 
be outlined which shows no particular complications in the development of a perturbative 
scheme around the hat null planes. 

First order vacuum theory can be obtained by setting the sources equal to zero, and 
taking u = 0. It is interesting to notice that under these conditions, fij satishes both the 
following: fij = 0 and fij = 0, from equations ([35|) and (^). This means that in vacuum 

our gauge conditions are a null version of the harmonic gauge. Furthermore, equation (p 


reduces to the homogenous wave equation DA = 0. Since all the components of fij can be 
found by applying partial derivatives and 9 or 8 operations to A (all of which commute with 
□) we have the result that every component of the metric satishes the homogeneous wave 
equation, Ofij = 0. We are currently exploring the possibility of hnding the general solution 
to (*), (A) and (C) in vacuum directly for A without the use of Z. This approach, which 
reverses the roles of Z and A as primary versus auxiliary quantities, and which yields new 
insights into the metricity conditions, will be presented elsewhere. 
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APPENDIX A: THE SOLUTION OF THE CUT FUNCTION EQUATION BY 

MEANS OF GREEN’S FUNCTIONS 


The equation for the linearized cut function (^) is of the form: 

m^z = s{x\cx) (Ai) 

where Z is a function of (C, Q and the four parameters S' is a given source. The solution 
to this equation is of the form: 

Z(C.C>“) = /G(C.C.C'.mc'.CVS' (A2) 


where G(C, C, C^ CO is the appropriate Green’s function for the operator 0^5^ on the sphere. 
This function, which has been found |p , is a rather simple expression in terms of the vector 
/“ introduced earlier in this work (see (|T^: 

1 


G(C,C,C',?) = j^i“(C,0'a(C',C')'n('‘(C,04(C',C')). 


(A3) 


It is also useful to have the explicit version in terms of (C, C); which is found from: 


^“(C,C)^a(C',C0 


(C-C0(C-C0 

(1 + CC)(1 + C'C0' 


There is a method for explicitly evaluating integrals of this type on the sphere for a large 
class of sources S p. This method consists of rewriting the integral as a sum of contour 
integrals around simple poles. In this way the integral is reduced in a large part to the 
calculation of residues. The method works for singular integrands as well, our present case. 
Explicit results have already been obtained by the authors for the vacuum case with arbitrary 
quadrupole radiation data of the form ct(m, C, C) = o('*^)2b^2m(C) 0) where 2b2m(C) C) is a spin- 
weight-2 spherical harmonic and a{u) is arbitrary. These results will appear elsewhere. 


APPENDIX B: CALCULATION OF THE DIVERGENCE AND DOUBLE 

DIVERGENCE OF THE METRIC 


For the component i 


only the use of (|2^) is needed: 


a>f-i = 3"/-o + s-/- 

= S“(i8Ai) + a-(-Ai) 

= ai(t3Ai)-a+(-Ai) 

= a,(t8Ai + An 

= 2^101/. 


(Bl) 


The component i 


0 requires the use of (^q) and the commutation relations: 


19 






^/o, =5%o + 5-/o- + a+/o+ 

= (9°(—-8^Ai) + d~fo- + diQdu + complex conjugate. 

= (9°(—-8^Ai) + (9“(-8 Ai) + (9i88z/ + complex conjugate. 

= (—(9°8 + (9“)-8Ai + (9i88z/ + complex conjugate. 

= —(i9i8 + (9+)-8 Ai + (9i88z/ + complex conjugate. 

= —(i9i8 + (9+)(2Sz/ — A+) + (9i88z/ + complex conjugate. 

= —2((9i 8 + (9+)8z/ + (9+((9i8 + (9+)A + diSBu + complex conjugate. 

= —2((9i 8 + (9+)8z/ + (9+(8(9i + 2d+)A + diBdu + complex conjugate. 

= —2((9i 8 + (9+)8z/ + (9+(4Sz/) + diBBu + complex conjugate. 

= —2(9i88z/ + 2(9+8z/ + diBBo + complex conjugate. 

= -2(9i88z/ + 2d+Biy + 2(9_8z/. (B2) 

The double divergence d'^d^ fij becomes: 

= d^d^foj + + d-d^f., 

= d,d^foj - d.d^u, - 

= (9i(—2(9i88z/ + 2d+Biy + 2d-Biy) — 2d-{diBiy) — 2d+{diBiy) 

= -2dlBBu. (B3) 


APPENDIX C: THE STRESS ENERGY TENSOR FROM THE NULL SURFACE 

POINT OF VIEW 

Since we have an explicit {(, () dependence of the component tu of a given symmetric 
trace-free tensor tij, the remaining components can be obtained by an appropriate combi¬ 
nation of successive applications of 8 and 8 to tu. The following have been obtained with 
the help of the vanishing trace condition, tu — toi + = 0; 

1 

ti- = -2^tn 
ti+ = -^Btu 
t— = ^BHu 
4+ = ^BHu 
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t+- = 

to- = - dtn 

to+ = - 8tii 

too = 


3 

^01 — ^ 99^11 + 2^11 


(Cl) 


As an example of the calculation, the hrst of ( |CI| ) was obtained in the following manner: 
til = Oid\tab = l^l’^tab therefore Stu = 2dl°'lHab = 2mf^ftab = —26°i6\tab = —2ti- and 
inverting: ti_ = — IStn. 

Since the components of the trace-free part of the stress tensor constitute the right side 
of the trace-free Einstein equations, then by the same operations suggested above for tij we 
can obtain the nine trace-free Einstein equations from dfu = tu. 


APPENDIX D: THE TRACE OF THE EINSTEIN SYSTEM AND THE 

METRICITY CONDITIONS 


Here we show that 


d(3niy + dlmiy) = 0 

by virtue of the metricity conditions. Though the calculation is fairly lengthy, we will outline 
the argument. We need Eq.(*), (A) and Eq.(|4^): 

□A = - 2did^u = -2{did - 2d_)du (Dl) 


(obtained from (*) and (A)). 

First we turn into a more suitable expression: 

ddiddu = did{did — 2d-)di> 9i(9+S -h 2di — 2do)Qi'. 


Now we use (pi|) to put the hrst term in terms of DA and use (*) to eliminate u in favor of 
A (also commute c?i9 in the hrst term): 

Qdlmu = -1 □ 8 Ai - i □ A+ + ^ai(a+g + 2^1 - 2ao)(A+ + i8Ai). 

If we put it together and after some work we obtain: 
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s(3 □ Z/ + = ^<9+(da + ((9i9 - 2(9_)(A+ + ^aAi)). 

If we commute B and d- in the term (9 _BAi and use (*) to eliminate A_ in favor of BAi we 
obtain: 

B(3 □ z/ + dlmv) = ^5+(da + (^iB - 29_)A+ + ^^iBBAi) 

—-(9iBBAi + -(9+9Ai — (9oAi + 2(9iAi^. 

With some more work using the commutation relations only we obtain: 

9(3 □ z/ + dlmv) = ^a+(-39+(A_ - ^9Ai) + ^(9i(9g - 99 + 4)Ai). 

This last result is clearly vanishing by virtne of the metricity condition (A) and the com- 
mntation relations for 9 and B on qnantities of spin weight eqnal to 2. 
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